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Abstract 

A theoretical and numerically study of dynamical properties in the sol-gel transition is pre- 
sented. In particular, the complex phenomenology observed experimentally and numerically in 
gelling systems is reproduced in the framework of percolation theory, under simple assumptions 
on the relaxation of single clusters. By neglecting the correlation between particles belonging to 
different clusters, the quantities of interest (such as the self Intermediate Scattering Function, the 
dynamical susceptibility, the Van-Hove function, and the non-Gaussian parameter) are written as 
superposition of those due to single clusters. Connection between these behaviours and the critical 
exponents of percolation are given. The theoretical predictions are checked in a model for perma- 
nent gels, where bonds between monomers are described by a FENE potential. The data obtained 
in the numerical simulations are in good agreement with the analytical predictions. 
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I. INTRODUCTION 



The gelation transition transforms a viscous liquid (sol) into an elastic disordered solid 
(gel). In general this process is due to the formation of a macroscopic molecule, due to 
the bonding of multifunctional monomers in solution, which makes the system able to bear 
stress. The extent of the gelation process may be measured by the monomer volume fraction 
<p, defined as = NV m /V, where N is the number of monomers, V m is the single monomer 
volume and V is the total system volume. On the static point of view, the sol-gel transition 
was interpreted 1,2 in terms of the appearance of a percolating cluster of monomers linked by 
bonds 3 , and experimental measurements of the geometric properties of gels have confirmed 
this correspondence (for a review see Stauffer et al. 4 and references therein). 

Complex dynamical behaviours are observed in gelling systems already in the sol phase. 
For example, light scattering measurements show non-exponential decay of the intermedi- 
ate scattering function, F(k,t), in both permanent 5 and thermoreversible physical gels 6 ' 7 . 
In particular power laws are observed at intermediate times, followed, at long times, by 
stretched exponential decays, exp(— (t/r) 13 ), with < < 1. In the gel phase, where ergod- 
icity is broken, only the power law decay survives. Usually the onset of stretched exponential 
decays (present also in other complex systems, as spin glasses and glassy systems) is associ- 
ated to the widening of relaxation times, which in gelling systems is due to the presence of a 
broad cluster size distribution close to the gelation threshold. However general predictions 
which connect this kind of relaxation to percolation theory are not easily feasible. 

In this paper, assuming "simple" behaviours for the relaxation of clusters with given size, 
we show how the "complex" phenomenology of the relaxation in permanent gels may be 
obtained from the superposition of the behaviours of clusters with different sizes. In partic- 
ular we are able to predict the behaviours of the self Intermediate Scattering Function, of 
the dynamical susceptibility, of the Van-Hove function, and of the non-Gaussian parameter, 
and to connect these behaviours to the cluster size distribution and the critical exponents 
of percolation. Then we check the theoretical predictions in a specific model for permanent 
gels, studied using Molecular Dynamics simulations. 

The paper is organized as follows. In Sect. II the analytical results are briefly summarized, 
and in Sect. Ill they are compared with the data obtained by Molecular Dynamics simulations 
of a model for permanent gels, where bonds between monomers are described by a FENE 
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potential 8 10 . In Sect. IV concluding remarks are discussed. Finally, in A, B and C the 
calculations are presented in details. 

II. CONNECTION BETWEEN STATIC AND DYNAMIC PROPERTIES 

In this section we summarize our calculations, which will be show in details in appendices. 
We consider a system of randomly distributed monomers with a fixed volume fraction, <fi. 
At time t = permanent bonds are introduced at random between monomers at a distance 
r < R, where R is suitably chosen. For a particular model see the FENE model 8-10 (Sect. 
IV), however the following arguments are independent on the details of the model. 

Following the percolation approach 1,2 , we identify the gel phase as the state where a per- 
colating cluster is present. We denote by <ft c the volume fraction of the percolation threshold. 
In our calculations we use some results from percolation theory 3 : In particular, in the sol 
phase, near the threshold, the cluster size distribution is given by n(s) oc s~ T e~ s / s * (where s* 
is a cut-off value given by £ D f , Df is the fractal dimension, and £ is the connectedness length 
which diverges at the threshold with the exponent u); in the gel phase, near the threshold, 
sn(s) is put equal to Pco5 s ,s max +Cs~ T+1 e~^ s ^ s *' ){d 1)/d , where is the density of particles in 
the percolating cluster of mass s max , d is the spatial dimension and C is a constant. More- 
over we assume that the relaxation time of clusters increases as a power law of the size 11 , 
r(s) oc s x . With these assumptions, in the hypothesis of simple behaviour for single cluster 
(exponential relaxation, simple diffusion, etc.), we obtain all the complex phenomenology 
observed experimentally and numerically near the threshold in gelling systems. In particular 
we are able to predict the behaviours of the self Intermediate Scattering Function, of the 
dynamical susceptibility, of the Van- Hove function, and of the non-Gaussian parameter. 

A. Self Intermediate Scattering Functions 

We first consider the self Intermediate Scattering Functions (ISF) : 

J^/(M) = KWM)>] (i) 
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where (...) is the thermal average over a fixed bond configuration, [. . . ] is the average over 
independent bond configurations of the system, 

1 N - 

** f (k,t) = ±Y,J k ' m) ~* m > ( 2 ) 

i=l 

and N is the number of particles. In the following we fix the wave vector k = k m i n and 
kmin = 27r/L, with L the linear system dimension. 

In terms of the contributions due to different clusters, F se if(k m i n ,t) can be written as 



Fself{kmim t) 



^ ^ stt(s) f s (k m i n , t) 



(3) 



where n(s) is the cluster size distribution (Nn(s) gives the number of clusters of size s) 



and f.(h,t) = (/c.(M)), where / C .(M) = J £ ieCs e ik <^~^ is the self ISF limited to a 
given cluster C s of size s, and ttt is the average over all clusters of given size s. 

By replacing the sum with the integral, the self ISF for a given bond configuration 
becomes: 

F^W)~/*»*)MW). (4) 

By assuming 

fs(kminit) — e : / T s ; (5) 

the integral in Eq.(4) gives in the thermodynamic limit the following predictions for the time 
dependence of F se if(k min ,t) in a permanent gel: 

(i) At the gelation threshold (0 = C ) 

F S elf{k mi n,t) OC t~ z T(z), (6) 

where = f da o z ~ x exp(— a) is the T-function with z = (r — 2)/x. 

(ii) In the sol phase (<fi < 4> c ) 

F se i f (k mmi t) cx r Cl e-^^, (7) 

where /3 = + 1), c\ = /3(r — 3/2), and r a oc (0 C — 0) , and / = xDfV. This 
approximated form, obtained in the long time limit, coincides with that suggested by 
Ogielski 13 as fitting function for the time dependent order parameter in spin glasses, 
and it is in agreement with experimental 5 and numerical 14 findings in gelling systems. 



(iii) In the gel phase (0 > C ) 

F se if(k mm ,t) ~ Poo + At'^e-^ , (8) 

where (3 = l/(x + l), and C\ = /3(r — 3/2) are the same exponents obtained in the sol phase. 
The plateau value, P^, gives the density of localized particles 15 . Clearly the main contribu- 
tion comes from localized particles of the percolating cluster, however a small contribution 
may be due to particles trapped inside it. 
The calculations are shown in details in A. 

Our findings given in Eq.s (6), (7) and Eq.(8) are in agreement with the theoretical 
predictions obtained in Ref. 16 in the Rouse and Zimm models for randomly cross-linked 
monomers, where x = 1 and x — 1/2 respectively. Similar calculations are also done in 
Ref. 17 in a different context. 



B. Fluctuations of the self ISF 



In Ref. 10 it was studied the dynamical susceptibility, defined as the fluctuations of the 
self ISF: 

X4 (M) = N [(\$ self (k,t)\ 2 ) - (WM)) 2 ] • (9) 

In particular it was shown that, in the sol phase, in the limit of t — > oo and k — > 0, X4,{k, t) 
coincides with the mean cluster size, S = s2n i s ) : 



lim lim Xi(k, t) = S, 

k— >0 t^rOO 



(10) 



which diverges at the threshold 10 with the exponent 7. Here we are interested in the time 
dependence of the dynamical susceptibility approaching the asymptotic value. 

We neglect the contributions due to disconnected particles at each time t. In this way 
we can write X4,{k,t) as a superposition of the contributions due to different clusters: 



x 4 (M) 



M*K/c.(M)#.(M)> - </<7.(M)X/S.(M)> 



(ii) 



where again fc s (k,t) = 1 J2 i€Ca e lk '^ 1 ^ r ^ ^ is the self ISF limited to a given cluster C s of 
size s, ttt is the average over all clusters of given size s, (. . .) is the thermal average over 
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a fixed bond configuration, and [. . . ] is the average over independent bond configurations. 
The term (fc s {k,t)fc (k,t)) in Eq.(ll) can be written as 

(fc s (k,t)f* Cs (k,t)) =1 £ ^Hm-rm^HrM-rM)) 

= ~2 Yl (e SKij{t) ), (12) 

where we have put Ay(t) = (ri(t) — fj(t)) — (r*j(0) — r}(0)). For connected particles i and 
j, \Aij(t)\ is finite, and in the low wave vector limit where |A,j(t)| -C 1/k we can assume 
e ik-A tj {t) ^ ^ Then, by supposing that (fc s (k,t))(f^ s (k,t)) = (fc s (k,t)) , in the zero wave 
vector limit the dynamical susceptibility for a given bond configuration can be written as: 

hmv 4 (M) * limJT^nOO (l - (f s (k,t)) 2 ) . (13) 

fe— >0 fc^O ' * 

s 

From this equation, using Eq.(5), it is direct to see that the dynamical susceptibility goes 
from zero (for t = 0) to the mean cluster size S (in the t — ■> oo limit), since the self ISF of 
clusters of given size, in the sol phase, goes from 1 (for t — 0) to zero (in the t — > oo limit). 

As in the previous section we can evaluate lim^o t) in the sol phase, < C . We 
find that, for time long enough, x^(k,t) approaches the asymptotic value in the following 
way: 

limx 4 (M) - S- (l - At C2 e- {2t/Ta)0 ^ , (14) 

where (3 = l/(x + 1) and c<i = /3(5/2 — r). The exponent (3 is exactly the same which 
appears in Eq.(7) for the decay to zero of the self ISF, the relaxation time in the stretched 
exponential function is given by r«/2, and finally the power law has a positive exponent C2 
different from the exponent C\ which appears in Eq.(7). 
The calculations are shown in details in B. 



C. Self part of the Van-Hove function 



The self part of the Van-Hove function is given by: 



G se if(r,t) = — 



N 



5>(r-|ri(0-ri(0)|)> 



i=l 



(15) 



If the motion of particles is diffusive with a diffusion coefficient D, G se if(r,t) = 
(l/4:iiDt) 3 / 2 exp(— r 2 /ADt), where r is the distance traveled by a particle in a time t. 
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Deviations from the Gaussian distribution were observed in different glassy and gelling 
systems 19,20 . In fact the van-Hove function seems fitted by a Gaussian only for short dis- 
tances, instead, for long distances, it is well fitted by an exponential tail that extends to 
larger distances for increasing times. 

The deviation from the Gaussian distribution indicates that some particles move faster 
than others, due to the presence of heterogeneities. In permanent gels, heterogeneities 
coincide with clusters of particles connected by bonds 10 . As matter of fact particles belonging 
to different clusters have a different diffusion coefficient depending on the cluster size. As a 
consequence it has been suggested 20 that, in the sol phase and in the diffusive regime (i.e. 
in the long time limit), G se if(r,t) is given by a superposition of Gaussians 



G se i f (r,t) 



sn(s 



(4vr J D(s)t) 3 / 2 



exp 



AD(s)t 



(16) 



where D(s) is the diffusion coefficient of clusters of size s and n(s) is the cluster size distri- 
bution. 

By assuming t _1 (s) oc D(s) = as~ x , and by replacing the sum with the integral in 
Eq.(16), predictions can be given for the dependence of G se if(r,t) on r and t. We find, in 
the limit r 2 ^> at: 

„2 



t 3 / 2 G self (r,t) oc 



1 xr 2 \ 1 

- 4 + F + 4^J exp 



r 
Aal 



(17) 



where A = 1 — r + 3x/2. 

The calculations are shown in details in C. 



D. Non-Gaussian parameter 

The non-Gaussian parameter is defined as 21 : 

3(Ar 4 (t)) 



a 2 (t) 



5«Ar 2 (t)»^ 



(18) 



where <Ar 2 (t)> = ^^x<]r,(t) -r,(0)| 2 >, and <Ar 4 (t)> = ± ElM(t) - r 4 (0)| 4 ). It is 
easy to show 21 that a 2 (t) is zero if the probability distribution of the particle displacements 



is Gaussian. 
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Using Eq.(16), in permanent gels the non-Gaussian parameter is expected to tend in the 
long time limit to a plateau, whose value is given by 



„ as 

a 





~D 2 -D 2 ~ 


- 1 = 






D 2 



(19) 



where, for each bond configuration, ttt is the average over the cluster distribution. From 
this relation it appears clear that the deviation from gaussianity («2 7^ 0) is due to the 
fluctuations of the diffusion coefficient, which in turns is related to the presence of dynamical 
heterogeneities, i.e. groups of particles with different diffusion coefficient. 

III. FENE MODEL FOR PERMANENT GELS 

In this section we check the theoretical predictions obtained in Sect. II in the FENE 
model for permanent gels. We consider a d — 3 system of N particles interacting with a soft 
potential given by Weeks- Chandler- Andersen (WCA) potential 22 : 

jjWCA = [ M(°/^) 12 ~ {o/r l3 f + |], r y < 2V*o 
13 \ 0, r tj > 2 1 /V 

where is the distance between the particles % and j. 

After the equilibration, at a given time t = particles distant less than i?o are perma- 
nently linked by adding an attractive potential: 

jjFene = I -0.5/c i?ohi[l - (ry/Ro) 2 ], < R 

[ OO, Tij > Ro 

representing a finitely extendable nonlinear elastic (FENE). The FENE potential was firstly 
introduced by Warner 8 and is widely used to study linear polymers 9 . We choose ko = 30e/er 2 
and Ro = 1.5a as usual 9 in order to avoid any bond crossing and to use an integration time 
step At not too small. 

We have performed Molecular Dynamics simulations of this model 10 : The equations of 
motion were solved in the canonical ensemble (with a Nose-Hoover thermostat) using the 
velocity-Verlet algorithm 23 with a time step At = 0.001<5r, where 5r = o{mje) x l 2 is the 
standard unit time for a Lennard- Jones fluid and m is the mass of particle. We use reduced 
units where the unit length is a, the unit energy is e and the Boltzmann constant ks is 
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FIG. 1: (Color online) Main frame: The self ISF, (f 3 {k,t)) for clusters of size s = 1,2,4, 10,21,52, 
for (f) = 0.09, k = k m i n . The curves are exponential functions, e */ T ( s ). Inset: The relaxation time, 
t(s), as function of the cluster size s for k = k m i n and <p = 0.09. The continuous curve is a power 
law s x with x ~ 0.65. 

set equal to 1. We choose periodic boundary conditions, and average all the investigated 
quantities over 32 independent configurations of the system. The temperature is fixed at 
T = 2 and the volume fraction <ft = 7ra 3 N/6L 3 (where L is the linear size of the simulation 
box in units of a) is varied from <fi = 0.06 to = 0.12. 

Using the percolation approach, we identify the gel phase as a state where a percolating 
cluster is present 1 ' 2 . With a finite size scaling analysis 10 we obtain that this transition is in 
the universality class of random percolation. In particular, we obtain that the cluster size 
distribution, n(s) ~ s~ T at the gelation threshold <p c = 0.09 ± 0.01, with r = 2.1 ± 0.2; the 
mean cluster size S((f>) = ^2 s 2 n(s) ~ (0 C — 0)~ 7 , with 7 = 1.8 ± 0.1; the connectedness 
length £ ~ (<f) c — cf))~~ v ', with v = 0.88 ± 0.01; and the fractal dimension of large clusters 
is Df = 2.4 ±0.1. In the following we fix the number of particles, N = 1000, where the 
threshold is <p c ~ 0.095. 
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FIG. 2: (Color online) Main frame: The mean squared displacement, (Ar 2 (s,i)}, for <j) — 0.09, 
and clusters of size s = 1,2, 4, 10, 21, 52. Inset: The diffusion coefficient, D(s), as function of the 
cluster size s for <f> = 0.09. The continuous curve is a power law s~ x with x ~ 0.67. 

A. Size dependence of dynamical behaviour of the clusters 

In the sol phase we have studied the dynamical behaviour of the clusters as a func- 
tion of the size s. In particular we have measured the self ISF and the mean squared 
displacement of clusters, respectively (f s (k,t)) = -J^ieC (e lk •( f K t ) -f i(°))) ) and (Ar 2 (s,t)) = 
sX^iec (l^(^) ~~ "^(O)l 2 )' where again ttt is the average of all clusters C s with fixed size s. 

After an initial transient, we find that {f s {k, t)) for k = k min is well fitted by exponential 
tail, e - */ 1 "^ (Fig.l), with 1/t(s) oc s~ x (Inset of Fig.l) and x not depending on the volume 
fraction. Our data furnishes 12 x ~ 0.65. 

The mean squared displacement of clusters (shown in the main frame of Fig. 2), after 
a ballistic regime at short time, displays a diffusive behaviour. The diffusion coefficient of 
clusters, D(s), obtained as (Ar 2 (s, t)) ~ 6D(s)t, is plotted in the inset of Fig. 2 as a function 
of size s. In agreement with the results for the relaxation time, we find that D(s) oc s~ x , 
with x ~ 0.67. 
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B. Self ISF and its fluctuations 



In the sol phase, due to the superposition of the contributions of different clusters, the 
self ISF is expected to follow Eq.(7): 



where /3 = 1/0+1) ~ 0.60, c x = /3(r-3/2) ~ 0.36, and r Q oc (0 c -</>) _/ , with / = xD f u ~ 
1.4. F se if(k min , t) is plotted in Fig. 3 for < <f> c . After the initial transient, the data are well 
fitted by the function, Eq.(7) (continuous curves in figures). Furthermore, in agreement with 
theoretical predictions, the relaxation time r a (plotted in the Inset of Fig. 3 as a function of 
0) appears to diverge approaching the transition threshold with the exponent / ~ 1.4. At 
the threshold <p c ~ 0.095 finite size effects appear. 

Interestingly the presence of dynamical heterogeneities, i.e. groups of particles with 
different diffusion coefficient, is related to the breakdown of Stokes- Einstein relation. As 
shown in A, the relaxation time r a is essentially the relaxation time of the critical cluster. 
On the contrary the diffusion coefficient of the system, obtained as (Ar 2 (t)) ~ 6Dt, is given 
by the average over clusters with different sizes, D = D = Y2 S sn(s)D(s). Since D(s) oc s~ x , 
it is clear that D is dominated by small clusters. As a consequence, although r a diverges at 
the threshold, D does not go to zero at <p c (see Inset of Fig.3), due to the diffusion of small 
clusters through the gel matrix even for <fi > <p c . 

In the gel phase a detailed analysis is not possible: F se if(k min ,t) (plotted in Fig. 4 for 
4> > (j) c ) displays a plateau, however at long times, due to finite size effects, it relaxes to zero 
following an exponential function. 

Finally in the sol phase, we have also measured the dynamical susceptibility, X4{k,t), 
defined as the fluctuations of the self ISF, given by Eq.(9). In Fig.5, X4(fc, t) is plotted for 
k = k m i n and different volume fractions. After the initial transient, the approach to the 
plateau is well fitted by Eq.(14): 



where (3 = l/(x + 1) ~ 0.6 and c 2 = /3(5/2 — r) ~ 0.24. Note that the relaxation time r' a 
coincides with r Q /2 only at low volume fraction; near the threshold r' a is instead lower than 
T a /2, due to the contribution of disconnected particles at intermediate times. 



Fselfikmimt) OC t 6 ^ ^ ^ , 



(22) 




(23) 
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FIG. 3: (Color online) Main Frame: Self ISF, F self {k,t) for k = k min and <j> = 0.06, 0.07, 
0.08, 0.085, 0.09, 0.095 (from left to right) as a function of time t. The lines are fitting curves: 
^-o.36 g -(t/r Q )°' 6 _ j nse t : Structural relaxation time, r a (full circles), compared with the inverse of 
the diffusion coefficient 20 /D (full triangles) as a function of the volume fraction. The full line is 
the fitting curve: t{4>) oc (4> c — <j))~ i with / ~ 1.4 and 4> c ~ 0.095. 



C. Self part of the Van-Hove function and the non-Gaussian parameter 

In the sol phase we have also measured the self part of the Van-Hove function, G se if(r, t), 
denned by Eq.(15). In the long time regime, G se if(r,t) is fitted by a Gaussian curve only 
for short distances, and it seems well fitted by an exponential function for long distances 20 . 
In the long time regime, clusters of any size show a diffusive behaviour (see Fig. 2), then 
we have suggested 20 that G se if(r,t) is given by a superposition of Gaussians, Eq.(16): 



G se i f (r,t) 



v— «, sn(s) r 2 



^ (ArrD(s)t) 3 / 2 6XP V 4D(s)t 



s 



(24) 



where D(s) is the diffusion coefficient of clusters of size s, and n(s) is the cluster size 
distribution. Data well agree with our hypothesis, as indicated in Fig.6, where we have 
used n(s) and D(s) measured in the simulations. As shown in the main frame of Fig. 7, 
e if(r, t) plotted as a function of r 2 /t for fixed <f> and for different times after the initial 
transient, collapse onto a single master curve, supporting the hypothesis that the data satisfy 
Eq.(24). 
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FIG. 4: (Color online) Self ISF, F self (k,t) for k = k min and <j> = 0.1,0.105,0.11,0.12. 




FIG. 5: (Color online) Dynamical susceptibility, XiiKnim t) f° r 4> = 0.06, 0.07, 0.08, 0.085 (from 
bottom to top) as a function of time t. The lines are fitting curves: S ■ ( 1 — j4t 0,24 e - ( i//r «)° 6 ) . 



Finally the comparison with the approximate form obtained in C for r 2 3> t, Eq. (17), 



t 3/2 G se// (r,t) oc 



xr 



r 1 


«2 


j ex p 


4at_ 



(25) 



(with A = 1— r+3x/2, s* and a obtained from the simulations) gives again a good agreement 
for high enough r (Fig. 7 and Inset). Note that the curves obtained from Eq.(25) and shown in 
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FIG. 6: (Color online) The self part of the Van-Hove distribution for <f> = 0.07 and time 
t= 93.199, 1285.02, 17715.2 (from top to bottom). Full lines are obtained from Eq.( 16). 

figure are numerically indistinguishable from exponential functions in the considered range 24 
(Inset of Fig. 7). 

In agreement with the above picture, we expect that the non-Gaussian parameter tends in 
the long time limit to a plateau, whose value coincides with the fluctuations of the diffusion 
coefficient, Eq.(19). In the main frame of Fig.8 the non-Gaussian parameter is plotted for 
different volume fractions, and in the inset of Fig.8 the plateau value is compared with 
the fluctuations of the diffusion coefficient. The data are in good agreement with Eq.(19), 
confirming that in permanent gels the non-gaussianity of the displacement distribution is due 
to the superposition of the contributions of clusters of different sizes. It is worth to notice 
that the main contribution to 0:2 (£) comes from small finite clusters. In fact, the bigger the 
cluster, the lower its diffusion coefficient and hence its contribution to the non-Gaussian 
parameter. Therefore, no criticality of the plateau value of the non-Gaussian parameter is 
observed approaching the transition threshold. 

IV. CONCLUSIONS 

In this paper we show how the complex dynamics, such as stretched exponentials and 
power law behaviors, observed experimentally and numerically in gelling systems, emerges 



14 




FIG. 7: (Color online) Main frame: t 3 / 2 G se// (r, t) as a function of r 2 /t for (f> = 0.07 and 
t = 93.199, 1285.02, 17715.2. The full line is obtained from Eq.(25). Inset: The self part 
of the Van-Hove distribution for 4> = 0.07 and time t = 93.199, 1285.02, 17715.2 (from top to 
bottom). Full lines are obtained from Eq.(25) with the values of s*, x, r and a measured in the sim- 
ulations. Note that the curves shown in figure are numerically indistinguishable from exponential 
functions. 
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FIG. 8: (Color online) Main frame: Non-Gaussian parameter, ct2(t), as a function of time t for 
<j) = 0.06, 0.07, 0.08, 0.09 (from bottom to top). Inset Asymptotic value of ot2(t) (empty squares) 
compared with the fluctuations of the diffusion coefficient given by Eq.(19) (full circles). 
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from the contribution of single clusters, which instead decay with a simple exponential. Fur- 
thermore, we establish a connection between this complex behaviour and critical exponents 
of percolation theory. 

We also find, in the diffusive regime, an asymptotic form (for long enough distances) 
of the self part of the Van-Hove function, which deviates from the Gaussian distribution, 
and is numerically very similar to an exponential tail, usually observed in a large variety of 
complex system. Our finding suggest therefore that such deviation from Gaussianity may 
be due to a general mechanism, which may be ascribed to the presence of heterogeneities. 

The theoretical predictions are found in agreement with numerical results, that we find 
in the FENE model for permanent gels. We suggest that a similar analysis can be extended 
to systems with finite lifetime bonds, as colloidal gels, glassy systems or spin glasses, where 
a "suitable" definition of clusters is necessary. 
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APPENDIX A: SELF INTERMEDIATE SCATTERING FUNCTIONS 



In this section we show in details the calculations which gives the predictions shown in 
Sect. Ill for the time dependence of the self ISF in the thermodynamic limit. 

Starting from Eq.(4), in the hypothesis that f s (k min ,t) ~ e - */ 1 "^, and 1/t(s) ~ as~ x , the 
self ISF for a given bond configuration becomes: 



F se i f (k min ,t) ~ / ds sn(s)e 



-ats 



(Al) 



Let us consider three different cases: fi] 



,11) (p < (p c ; (111) <p > 



For <f) < (p c , n(s) can be written as n(s) ~ s T exp(— s/s*) 3 , where s* is a cutoff value 25 . 
Then we obtain: 



F se i f (k min ,t) ^ J ds exp(t(s)), 



(A2) 



where 



t(s) = -fr - l)lns - — - — . 



(A3) 

s" s* 

(i) At the gelation threshold, = <f) c , s* — > 00, and the integral, Eq.(A2), with t(s) given 
by Eq.(A3), can be evaluated exactly: 



Fself(k m in,t) OC t Z T(z), 



(A4) 



where T(z) = J °° da a z ~ x exp (—a) is the T-function with z — (r — 2)/x. 

(ii) In the sol phase, <fi < <fr c , we are able to give only approximated predictions. The 
function t(s), given by Eq.(A3), has a maximum for s such that 26 



- = -0—1 

s 



xat 

s x 



Let us approximate t(s) with t(s) — (s — s) 2 /(2a 2 ), where 



d 2 t(s) 







ds 2 



^(-(t-1)+x(x + 1)^\ 



If s ^> a, J ds exp 



(fzg)! 

2<7 2 



(2tt ( t 2 ) 1 / 2 , and 



oc - 



(t-2) 



F se if{k min ,t) ~ exp [t(s)] / (is exp 



[r-l) + x(x + l) — 



fs - s) 2 



2a s 





/ s at\ 


1 exp 




~ \s* + F/_ 



(A5) 



(A6) 



(A7) 
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F self (k min ,t) <x g * J ( x(t - 1) + (x+ 1)- ) (A8) 



Let us consider two limit cases: (1) s* — > oo; (2) at/s x —* oo. 

(1) Using Eq.(A5), Eq.(A7) can be written in the following way: 

exp |^(^-i) + (*+i)| w/2 

which in the limit s* — » oo, where 5 ~ (xat/(r — 1)) , gives again F se if(k min ,t) oc t 
with z = (r — 2) /rr, in agreement with previous calculations. 

(2) Using Eq.(A5), Eq.(A7) can be written in the following way: 

(x+l)at 



cxp 



s 



a^^- 1/2 



F se i f (k min ,t) oc ^—^ J - [ -(t-1) + x(x + 1)— ) (A9) 



In the limit at/s x — > oo, s ~ (xas*tY^ x+1 \ and we obtain 

F self (k min ,t) oc r Cl e- (i/Ta) ", (A10) 

where (3 — l/(x + 1), c\ — (3{r — 3/2), and r a oc s* x , which diverges at the threshold as 
power law with the exponent / = xDju. 

(hi) Finally in the gel phase, <fi > <fr c , sn(s) in Eq.(Al) can be written as 3 ' 27 sn(s) = 
Poo& 's, Smax +C 's~ T+1 e~( s > ' s *^ d d , where is fraction of particles belonging to the percolating 
cluster, s max is infinite in the thermodynamic limit, d is the spatial dimension, and C is a 
constant. Then, from Eq. (Al), we obtain: 

F se i f (k min ,t) -Poo + C J ds exp(t(s)), (All) 

where 

/on (d-l)/d nf 

t(s) = -(T-l)lns- (±Y (A12) 



,S / s 

Following the same arguments as in the sol phase, the second term in Eq. (All) is written 

as: 

xd \ / at\ 



J ds exp(t(s)) oc s ( T 2 -* exp 



d - 1 J \s x ) 



,,f I+1 ')^, (A13) 



d \ d J s 

where the maximum point, s, of t(s), Eq. (A12), is obtained from 



axt d — 1 / s 



(d-l)/d 

+ -') + 7 = Ti?l (A14) 
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In the limit at/s x — > oo, s oc tVOH- 1 and we obtain 

Fselfi^mini ty — Poo ~t~ At ® C ^ ^ 



(A15) 



where /3 g = (1 - l/d)/(l - 1/d + x), = /3 fl (r - 3/2 - l/2d)/(l - 1/d). 

For a finite system however s max is finite, and the behavior of F se if(k m i n ,t) is given by: 



- i oo c - v ^ i max j 



where r, 



mas: ^"l^maa; 



Fself {feminity — 

) is the relaxation time of the percolating cluster. 



(A16) 



APPENDIX B: FLUCTUATIONS OF THE SELF ISF 



In the hypothesis of the previous section (i.e. (f s {fe,ty) — e */ r ( s ) and l/r{s) ~ as x ) 
Eq.(13) becomes 

lim X 4(A;,t) = ^s 2 n(s)(l-e- 2t/r(s) ), (Bl) 



which, for < C , by replacing the sum with the integral, gives: 



lim x±{fe, ty — S — I dsexp{w{s)), 

k — >0 



where 

s 2at 

wis) = — (t — 2) In s . 

v i \ i s * s x 

Note that the function w{s) has a maximum for s such that 

— = - (t - 2 + . 

s* v ; s x 

Let us approximate u>(s) with w{s) — (s — s) 2 /(2cx 2 ), where 



a 



1_ = _ d 2 w(s) 

~2 = 



ds 5 



=2 |(r-2)-x(x + l) 



If s ^> a, J (is exp 



(2 7 ra 2 ) 1 / 2 , and 



(s - s) 2 

5 — lim x 4 (/c, t) ~ exp [w(s)] / dsexp ( — 



r- 2) -*(* + !) — J 



Y 1 ' 2 




1 exp 






\s* s / 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 
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The limit at/s x — > oo, where s oc (xs*ty^ x+1 \ gives 

lim Xi(k, t) ~ 5 • (1 - At C2 e- {2t/Ta ^), (B7) 

where /3 = l/(x + 1) and C2 = /3(5/2 — r). The exponent /3 is exactly the same which appears 
in Eq.(7) for the decay to zero of the self ISF. Note that the relaxation time in the stretched 
exponential function is given by r a /2. 



APPENDIX C: SELF PART OF THE VAN-HOVE FUNCTION 



In the sol phase, where after an initial transient the system is found in a diffusive regime, 
we assume the validity of Eq.(16) for the self part of the Van-Hove function. By replacing 
the sum with the integral, and putting D(s) ~ as~ x , G se if(r,t) can be written as 

1 



G self (r,t) 



(4vrt) 3 / 2 



ds exp (z(s)) 



where 



z s 



1 - T 



3x\ s 
- In s 



s x r 2 



2 J s* Aat 

The condition for the first derivative of z(s) to be zero gives: 



2 T 

s* Aat 



1-T+-X. 



(CI) 



(C2) 



(C3) 



This equation admits a solution s m > only if A = 1 — r + 3x/2 > 0. Under this 
condition the solution of Eq.(C3) is a maximum, in fact 



d 2 



zis 



ds 2 



— [A + x(x - 1) 



2 n x 



r s 



Aat 



>0, 



(C4) 



is always true for A > (i.e. x > 2(r — l)/3) 28 . In this case we can approximate z(s) 



with z(s m ) - (s - s m ) 2 /(2cr 2 ), where ^ 
(27rcr 2 ) 1//2 , and we can write: 



d 2 z(s 



If s m ^ a, J ds exp 



2a l 



t 3/2 G self (r ,t) oc exp [z(s m )] / ds exp 



2a 2 



oc s 



l+A 



exp 



Aat 



A + x{x- 1) 



r2 S X m 

Aat 



-1/2 



(C5) 



Let us consider two limit cases: (i) r 2 <C at, and (ii) r 2 ^> at. From Eq.(C3) we obtain: 

As* r 2 < at 



' 4i ^ Vl r 2 »at, 



(C6) 



20 



and Eq.(C5) becomes: 



t^G self (r,t) oc { V + 

(4Aat\( ) /X 
\ xr 2 ) 



, -1/2 


r 2 " 


J exp 


4o't 



exp 



s* \ r 2 / 



r 2 < at 
r 2 > at, 



(C7) 



where a' = a/(As*) ;E , and A' = (A/x)^^ x \ Note that the condition s m 3> cr is satisfied in 
the limit r 2 /at — > where s m — > As*, which increases with increasing x and/or approaching 
the gelation threshold. In the opposite limit, r 2 /at — > oo, s m — > 0, and hence the condition 
s m > a is not satisfied. In this case this approximation is expected not to hold and a 
development at the first order of the Taylor series of z(s) around s — 1 should be more 
appropriate (see below). 

In the case A = 1 — r + 3x/2<0 (as in the FENE model for permanent gels presented 
in Sect. Ill), z(s) is a monotonic decreasing function of s: 

ds s s* Aat 

Let us develop z(s) at the first order of the Taylor series around the maximum s m = 1: 

dz 



z(s 



Z(S r 



+ 



ds 



IS - S r 



(C9) 



By replacing Eq.(C9) in Eq.(Cl), we obtain 



t 3/2 G self (r,t) oc 



dz 
ds 



exp [z(s m )] 



oc 



1 xr 2 \ 1 

A + 7 + ^) exp 



r 

Aat 



(CIO) 



This approximation is expected to hold in the limit of long distances (r 2 ^> at), where the 
second order of the Taylor series around s m = 1 is much smaller than the first one, and only 
a small interval of values of s around s m = 1 contribute to the integral in Eq.(Cl). 

The research is supported by CNR-INFM Parallel Computing Initiative, S.Co.P.E. and 
L.R. N.5 2005. 
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"m 

d 2 z{t 



28 By writing d ^ 



ds 2 



we easily see that, for A > 0, d 



1 / ly^ 1 o*^ 



,2 
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